Charge and spin correlations of a one dimensional electron gas on the continuum 
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(Dated:) 

We present a variational Monte Carlo study of a model one dimensional electron gas on the 
continuum, with long-range interaction (1/r decay). At low density the reduced dimensionality 
brings about pseudonodes of the many-body wavefunction, yielding non-ergodic behavior of naive 
Monte Carlo sampling, which affects the evaluation of pair correlations and the related structure 
factors. The problem is however easily solved and we are able to carefully analyze the structure 
factors obtained from an optimal trial function, finding good agreement with the exact predictions 
for a Luttinger-like hamiltonian with an interaction similar to the one used in the present study. 



I. INTRODUCTION 

Electrons in low dimensions, as found in modern semi- 
conductor devices are greatly affected by correlation 
effects that may dramatically change their behavior and 
bring about new phenomena and phases. Technological 
advances have reached the point where it is now possible 
to fabricate high mobility systems in which the electrons 
move essentially free and are intentionally confined at the 
nanoscale in suitable geometries. Though the details of 
the devices still play a role at the quantitative level0, to 
a good approximation some of these systems provide a 
very close realization of simple electron gas models. We 
should emphasize that this is true only in low dimension 
and the experimental realization of a clean 3-dimensional 
(3D) electron gas (EG) in the quantum regime is in fact 
still out of sight, as is the long sought 3D Wigner crystal. 
Here, we focus our attention on a simple ID model with 
long-range interaction, i.e., a quantum wire, with the aim 
of examining its structural properties. 

A large variety of models have been proposed over 
the years for quantum wires. One may roughly distin- 
guish two classes of systems, depending whether the wire 
is obtained from the confinement of a 3D|j, Q or of a 
21)0 00 electron gas; or one may focus on the nature of 
the confinement, which can be harmonic [JQ- hard- wall 
like 0, El or CoulombicQ- In all these cases, the effec- 
tive electron-electron interaction in the lowest subband of 
the transverse motion turns out to have a 1/r long-range 
tail, whereas the details of the confinement only affect 
the behavior at short range, which displays anyhow a di- 
vergence milder than 1/r, if any. Thus, if the interest is 
in the effect of the long range tail of the interaction, one 
may well resort to ad-hoc choices of the interaction 
or even to tight-binding models yielding Hubbard-like 
Hamiltonians with long range intersite interactional fiof : 
however, this last choice may introduce^] processes 
(such as the Umklapp scattering) not present in contin- 
uum models. One should also recall that, for interacting 
Fermions in ID, the familiar Fermi-liquid paradigm must 
be abandoned in favor of the Tomonaga-Luttiger liquid 
concept pd|. 

The above model wires have been studied at T = 



employing several techniques ranging from DFT 
treatments |3 to dielectric- formalism approaches of the 
STLS type|a 0, 

from exact diagonalization|l3l to 
density matrix renormalization group techniques 9] , to 
bosonization techniques@ _ the last approach being able 
to provide exact results for some properties if a lin- 
earized kinetic energy dispersion is assumed. Investiga- 
tions of the homogeneous 1DEG on the continuum, per- 
formed with either the variational (VMC) or the diffu- 
sion Monte Carlo (DMC), are also available [HE EU . 
Contrary to STLS predictions [T^ for the same model, 
DMC simulations [l4L Hq yield no Bloch instability (spon- 
taneous magnetization of the system) , in accordance with 
the Lieb and Mattis theorem[llj . They predict pair cor- 
relations that are consistent with a quasi Wigner crys- 
tal (strong charge correlations resulting in a divergent 
peak in the charge structure factor at 4kp 0). However, 
they also appear to imply the presence of quasi antifer- 
romagnetic order, i.e., a divergence of the spin structure 
factor at 2Uf in the unpolarized state, a result at vari- 
ance with the findings of the bosonization technique 0. 
Rather than to the differences between models [ljij], such 
a discrepancy is likely to be due to the ergodicity prob- 
lems encountered in both DMC and VMC simulations 
at strong coupling^, 0, as exchange between oppo- 
site spin electron is frozen out by the strong Coulomb 
repulsion. While in the DMC method a solution of 
such a problem is not straightforward 16], in VMC sim- 
ulations one may device fairly simple ways of restoring 
ergodicitv [hH Hq . thus establishing the true nature of 
spin correlations, as we show below. 



II. MODEL AND WAVEFUNCTION 

The Hamiltonian of our model quantum wire|3j is, in 
units in which <Zq = h 2 e/(m*e 2 ) = 1 and e 2 /(2eajj|) = 1, 

N 

i=l i<j 

V b (x) = (V^/26)exp (x 2 /4b 2 ) erfc {\x\/2b) , (2) 
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with b a measure of the wire width. Above, e is the 
dielectric constant of the semiconductor and m* the ef- 
fective mass of the carriers. The pair interaction Vb(x), 
which is finite at the origin, V&(0) = y / 7r/(26) and de- 
cays as 1/x for x 3> b, is obtained from an harmonic 
confinement of the 3DEG, after projection on the lowest 
subband of transverse motion|3j; it thus provides a good 
approximation to the 3D system at low density p, i.e., 
r, = l/2p ^> 7r6/4, with the density parameter r s also 
providing an estimate of the Coulomb coupling, as ratio 
between average potential and kinetic energies. 

We use VMC simulations and focus on the ground 
state properties of a thin wire (b — 0.1) on a fairly 
wide coupling range (1 < r s < 10), thus considering 
only unpolarized states[l5l To this end, we re- 

sort to a Slater-Jastrow wavef unct ion [l8j — JD^D^, 
with D" a determinant of N a plane waves and J = 
ex P[~ Si<i < u (\ x i ~ x j\)]i with u(x) a two-body pseu- 
dopotential to be optimized. This is the simplest cor- 
related wavefunction for an unpolarized Fermion state 
with homogeneous density and can be further improved 
with the inclusion in the Jastrow factor J of higher order 
pseudopotentials^^. Our goal here is a careful finite- 
size scaling extrapolation of structural properties, to ob- 
tain the long-range behavior of pair correlation functions. 
Thus, to reduce the finite size bias we keep N~= odd, 
in order to avoid shell degeneracy effects, use periodic 
boundary conditions and Ewald-sum |14l Ha. Il8| the pair 
potential Vb(x). 

We remind the reader that in ID the nodes of the 
Fermionic ground state (of given symmetry) are know 
exactly 20], as they are fully determined by exchange 
(anti)symmetry, and coincide in particular with those of 
the above wavefunction for the unpolarized state. Thus, 
DMC provides in principle exact energies [Til ITU IT^| that 
may serve as benchmark for the VMC simulations when 
optimizing the Jastrow factor J. As starting point for 
our simulations we take a two-body pseudopotential of 
the RPA type EE), which in Fourier space reads 

2pu RPA (k) = So(k)- 1 + /%(fc)- 2 +4pVUfc)/fc 2 , (3) 

with5 (fc) = (k/2k F )8(2k F -k)+0(k-2k F ) the structure 
factor of a non interacting 1DEG, 

V b {G)=E 1 (b 2 G 2 )e X p(b 2 G 2 ) (4) 

the Fourier transform of the interaction, and 9(x) and 
Ei(x) respectively the Heaviside and the exponential in- 
tegral functions. 

The repulsive nature of the pair interaction Vb(x) is 
directly reflected in the pseudopotential, which is shown 
in Fig. ^ The pseudopotential is repulsive and the re- 
pulsion increases appreciably with decreasing the den- 
sity (increasing the r s ). Thus, as the density is lowered 
electrons are kept apart more and more effectively and 
this results in a quasi long range order which can be de- 
scribed as a quasi Wigner Crystal 8]. Though the pseu- 
dopotential remains finite at contact, the Jastrow factor 
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FIG. 1: RPA pseudopotential for b = 0.1 and r s = 1 (solid 
line), r s = 2 (dashed line), r s = 6 (dotted line), r s — 10 
(dot-dashed line). 



becomes exponentially small yielding what may be seen 
as pseudonodes of the wavefunction. These pseudonodes 
have no particular effect on like spin electrons, as the 
wavefunction is already vanishing at contact for such elec- 
trons and most importantly particle exchanges are explic- 
itly summed over in the determinants. On the contrary, 
the effect on opposite spin electrons when combined with 
the reduced dimensionality may become dramatic. In 
a random walk in configuration space with importance 
sampling given by the Slater-Jastrow wavefunction, the 
RPA pseudopotential tends to freeze out the exchange 
between opposite spin electrons and may cause ergodic- 
ity problems, which show up in spin correlations |l4llT6| . 
Evidently, it is not only the presence of the pseudon- 
odes to cause problems but also the slope with which the 
pseudonodes are approached, on the scale of the inter- 
particle distance. When such a slope becomes sufficiently 
large, naive algorithm may become inefficient in sampling 
inequivalent pockets in configuration space, delimited by 
the pseudonodes. 



III. ERGODICITY AND VMC OPTIMIZATION 

Ergodicity becomes an issue in the present context be- 
cause in ID there is no exchange without crossing, for 
opposite spin electrons. With increasing the coupling, 
the maximum displacement 5 has to be made smaller (on 
the r s scale), to keep the acceptance high, and this re- 
duces the exchange frequency v = spin flips/total moves 
substantially. The problem is very serious in DMC simu- 
lations, as very small displacements are necessary |22j, to 
keep the time step error small, and a direct sampling of 
spin exchange is not straightforward. On the contrary, in 
the VMC method it is relatively easy to devise effective 
ways of keeping v above a suitable threshold (~ 0.01) 
ensuring ergodicity. The simplest choice is to work at 
acceptance appreciably lower than 50%, i.e., at larger S . 
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Sampling 
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-0.474825(9) 
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-0.474827(9) 


*G 


1.3 
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7.36 • 10 _;i 


-0.474832(8) 



TABLE I: Exchange frequency v at r s = 10 and TV = 22. 5, 
A and E are the maximum displacement of the MC move, the 
average acceptance, and the energy per particle. In all cases 
an optimized (scaled) RPA pseudopotential (see text) is used. 



Other possible choices involve either umbrella sampling 
/reweightingp^ . i.e. the use of configurations generated 
from a guidance function less repulsive than the trial 
function considered above, or the use of a general- 
ized Monte Carlo move, involving the direct sampling of 
exchange between opposite spin electrons |l4j. One may 
also consider combinations of the above. In Tab. we 
focus our attention on the first two choices, for the worst 
case examined here, i.e., r s = 10. It is clear that both the 
large <5 and the recipes are effective in keeping v close 
to value ss 0.01 found at small couplings, r a « 1, where 
there are no apparent signals of nonergodicity. Indeed we 
find that keeping v w 0.01 the noisy patterns previously 
found in the spin structure factor [14l Il6| disappear. On 
the other hand, the energy does not seem to depend on 
v and this is likely due to the fact that at strong cou- 
pling different spin configurations are almost degenerate. 
In the following we restrict our attention to the large S 
prescription which has the advantage of maximum sim- 
plicity. 

To get an accurate description of correlation functions 
we systematically optimize ^t, employing the variance 
minimization method |24j. As the plane-wave determi- 
nants provide the exact nodes for our ID system, we only 
need to optimize the Jastrow factor J, i.e., the two-body 
pseudopotential u(x), for which we have considered both 
a scaled RPA form u a {x) = otURpA{x) and a systematic 
expansions in terms of Chebyshev polynomials 

u{x) = J2 a m T 2m (2\x\/L - 1), (5) 

ra— 1 

satisfying the continuity of the first derivative at the edge 
of the simulation box, u'(L/2) = 0. The results of the 
optimization are summarized in Tab. \n\ in terms of the 
correlation energy E c = Eq — Ehf, with Eq and Ehf the 
exact (DMC) and Hartree-Fock energies. We find that 
the scaled RPA form is essentially equivalent to the true 
two-body pseudopotential, obtained with about 10 terms 
in the sum of Eq. (0. Moreover, adding a three-body 
term in J, within a factorization ansatz a la Feynman 
retaining the first two /-components docs not improve 
our results in any appreciable manner. Therefore in the 
following we use a scaled RPA psudopotential optimized 
for each N. 
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FIG. 2: Peak in the charge structure factor versus the system 
size. The curves are fits to the VMC predictions (symbols). 

IV. PAIR CORRELATIONS 

According to the predictions of the bosonization tech- 
nique applied to a Luttinger Hamiltonian (linearized ki- 
netic energy) with long-range interaction 8] , charge cor- 
relations exhibit long range tails resulting in a divergent 
peak at AUf in the charge structure factor S++(k) = 
(p + (k)p + (—k)) /N, whereas spin correlations decay faster 
and a finite peak is found at 2k p in the spin structure fac- 
tor S— (k) = (p-(k)p-(-k))/N. Here, p±{k) = Pt(k)± 
Pl(k) is the Fourier component of the charge (spin) den- 
sity. We have performed VMC simulations for b = 0.1, 
r a = 1, 2, 6, 10, N = 10, 22, 42, 82, 162, and computed the 
above structure factors at several k values. In Fig. Inland 
13 we give the peak heights N{Akp) = S ++ (4kp, N) and 
S(2kp) = S--(2kp, N), which appear to be increasing 
functions of N. As such heights are determined by the 
long-range behavior of correlations we would expect them 
to scale with N similarly to what found for the long-range 
Luttinger liquid of [8j, for which at large N one readily 
obtains 

S(2k F ) = a 3 ^\ogL + l/c) exp(-cVlogL) + a 4 , 
N(Akp) = ai£cxp(-4cVlog£) +a 2 , (6) 

with L — 2r s N and c, 01,02,03,04 model and density de- 
pendent. We have used Eqs. El to fit our results, obtain- 
ing first ai,a2,c from the fit of N(4kp) and then 03,04 





Etot 


E c 


%E C 


RPA 


-0.47207(2) 


-0.20519(56) 


0.9830(15) 


Scaled RPA 


-0.474825(9) 


-0.20794(55) 


0.9962(15) 


Chebyschev 


-0.474900(9) 


-0.20802(55) 


0.9965(15) 



TABLE II: Total energy Etot, correlation energy E c and per- 
centage of correlation energy %E C for b = 0.1, r s = 10, and 
N = 22. The fraction of the correlation energy recovered is 
computed from DMC calculations which provide the exact GS 
energy for a 1DEG. 
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FIG. 3: Peak in the spin structure factor versus the system 
size. The curves are fits to the VMC predictions (symbols). 
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= 1 


4.1 


3.5 


0.71 


r a 


= 2 


4.7 


1.5 


0.78 


Vs 


= 6 


0.55 


0.03 


0.80 


r s 


= 10 


3.6 


1.3 


0.88 



TABLE III: Fit of the peak of the charge and spin structure 
factors: reduced \ 2 an d interaction (fit) parameter c. 



from the fit of S(2k F ), with c fixed by the fit of N(4k F ). 
As it can be seen from Fig. [3 01 and Tab. II I II the long- 



range tails implied by Eqs. HJare indeed consistent with 
our VMC results, yielding a finite S(2k F ) in the thermo- 
dynamic limit. We have reason to believe that a further 
improvement of the \ 2 would be obtained by discard- 
ing the smaller values of N in favor of larger ones, on the 
ground that Eqs. Elonly embody the long-range behavior 
of correlations. 



V. SUMMARY 

In conclusion, we have shown that a relatively simple 
(unbiased) optimization of the wavefunction is able to 
yield correlations functions for a 1DEG on the contin- 
uum which are in qualitative agreement with the exact 
predictions of Schulz 8] for a Luttinger Hamiltonian with 
long-range interactions. The correlations that we find do 
support the existence of a quasi Wigner crystal and a 
non-ordered unpolarized magnetic phase. These results 
are at variance with what previously found for the same 
model wire with VMC HI and DMC0 simulations, 
as those calculations were biased by non-ergodicity. The 
results found here are strongly corroborated by very re- 
cent simulations 25| with a novel lattice regularized DMC 
(LRDMC)method|26| . which significantly alleviates the 
lack of ergodicity present in the standard DMC. In partic- 
ular we find 25] that the time projection implemented by 
the LRDMC only produces minor quantitative changes 
in the charge and spin peak heights, leaving our present 
conclusions unchanged. 



[1] Nanostructure Physics and Fabrication, ed. by M. Reed 
and W.P. Kirk (Academic Press, Boston, 1989); C. Weis- 
buch and B. Vinter, Quantum Semiconductor Structures, 
(Academic Press, Boston 1991). 

[2] See, e.g., S. De Palo, M. Botti, S. Moroni, and G. Sena- 
tore, cond-mat/0410145 and to be published. 

[3] W. I. Friesen and B. Bergesen, J. Phys. C 13, 6627 
(1980). 

[4] A. Gold and A. Ghazali, Phys. Rev. B 41, 7626 (1990). 
[5] G. Y. Hu and R. F. O'Connel, Phys. Rev. B 42, 1290 
(1990). 

[6] G. Y. Hu and S. Das Sarma, Phys. Rev. B 48, 5469 
(1993). 

[7] Y. Sun and G. Kirczenow, Phys. Rev. B 47, 4413 (1993). 
[8] H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 
[9] G. Fano and F. Ortolani, A. Parola, L. Ziosi, Phys. Rev. 
B 60, 15654 (1999). 
[10] S. Capponi and D. Poilblanc, T. Giamarchi, Phys. Rev. 

B 61, 13410 (2000). 
[11] K. Schonammer and V. Meden, Am J. Phys. 64, 1168 

(1996); J. Voit, Rep. Prog. Phys. 58 977 (1995). 
[12] L. Camels and A. Gold, Europhys. Lett. 39, 539 (1997). 
[13] D. Poilblanc, S. Yunoki, S. Maekawa, and E. Dagotto, 
Phys. Rev. B 56 R1645 (1997). 



[14] A. Malatesta, Ph.D. Thesis (in english), University of Tri- 
este (1999); it can be obtained from senatore@ts.infn.it. 

[15] A. Malatesta and G. Senatore, J. Phys. IV 10, Pr5-341 
(2000). 

[16] M. Casula, Laurea Thesis (in italian), University of Tri- 
este (2001); it can be obtained from senatore@ts.infn.it. 

[17] E. Lieb and D. Mattis, Phys. Rev. 125, 164 (1962). 

[18] See, e.g., D. Ceperley, Phys. Rev. B 18, 3126 (1978). 

[19] See, e.g., S. Moroni, S. Fantoni, and G. Senatore, Phys. 
rev. B 52, 13547 (1995). 

[20] D. M. Ceperley, J. Stat. Phys. 63, 1237 (1991). 

[21] T. Gaskell, Proc. Phys. Soc. 77, 1182 (1961), ibidem 80, 
1091 (1962). 

[22] P. J. Reynolds and D. M. Ceperley, B. J. Alder, W. A. 

Lester,Jr., J. Chem. Phys. 77, 5593 (1982). 
[23] D.M. Ceperley and M. Kalos, in Monte Carlo Methods in 

Statistical Physicsm ed. K. Binder (Springer, New York, 

1979). 

[24] C. J. Umrigar, K. G. Wilson, J. W. Wilkins, Phys. Rev. 

Lett. 60, 1719 (1988). 
[25] M. Casula, S. Sorella, and G. Senatore, to be published. 
[26] M. Casula, C. Filippi, and S. Sorella, |concTm at /0502388 

and to be published. 



